Tensor-optimized shell model for the Li isotopes 
with a bare nucleon-nucleon interaction 
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We study the Li isotopes systematically in terms of the tensor-optimized shell model (TOSM) 
by using a bare nucleon-nucleon interaction as the AV8' interaction. The short-range correlation is 
treated in the unitary correlation operator method (UCOM). Using the TOSM+UCOM approach, 
we investigate the role of the tensor force on each spectrum of the Li isotopes. It is found that 
the tensor force produces quite a characteristic effect on various states in each spectrum and those 
spectra are affected considerably by the tensor force. The energy difference between the spin-orbit 
partner, the p\/2 and P3/2 orbits of the last neutron, in 5 Li is caused by opposite roles of the tensor 
correlation. In 6 Li, the spin-triplet state in the LS coupling configuration is favored energetically 
by the tensor force in comparison with jj coupling shell model states. In 7,8,9 Li, the low-lying 
states containing extra neutrons in the p 3 / 2 orbit are favored energetically due to the large tensor 
contribution to allow the excitation from the Os-orbit to the P1/2 orbit by the tensor force. Those 
three nuclei show the jj coupling character in their ground states which is different from 6 Li. 

PACS numbers: 21.60.Cs, 21.10.-k, 27.10. +h 27.20. +n 



I. INTRODUCTION 



It is an important subject in nuclear physics to un- 
derstand the nuclear structure from the viewpoint of the 
nucleon-nucleon (TV TV) interaction. The NN interaction 
has strong tensor forces at long and intermediate dis- 
tances caused by the pion exchange and strong central 
repulsion at short distance caused by the quark dynam- 
ics. It is important to investigate the nuclear structure by 
treating these characteristics of the NN interaction [j]- 
0| ■ Recently, it became possible to calculate nuclei up to 
mass around A ~ 12 using a NN interaction with the 
Green's Function Monte Carlo method (GFMC) HQ. 
At present, this method requires the extreme computa- 
tional time to be applied to heavier nuclei. It is desired 
to develop a new method to calculate nuclear structure 
with large nucleon numbers by taking care of the char- 
acteristic features of the NN interaction. 

The presence of the tensor force in the NN interaction 
induces d-wave components in a nucleus, in particular, 
for proton-neutron (pn) pairs as the deuteron. The d- 
wave component in the deuteron makes the system bound 
via the large s-d coupling of the tensor force. This d- 
wave component is also found to be spatially compact as 
compared with the s-wavc component due to large mo- 
mentum characters brought by the tensor force |6[ . This 
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effect originates from the pseudo-scalar nature of the one- 
pion exchange. It is reported experimentally that a large 
fraction oipn pairs is observed as compared with pp or nn 
pairs in light nuclei 0, Q ■ This enhancement of pn pairs 
is hard to reproduce theoretically in a simple shell model 
(mean-field picture) Hexcept for the rigorous method 
as the one of GFMC [9|, which treats the tensor force 
explicitly. It is important to make efforts to study high 
momentum components caused by the tensor force in fi- 
nite nuclei (loj . 

There are two important developments to perform nu- 
clear structure calculations in heavy nuclei by including 
the necessary dynamics induced by the bare NN inter- 
action. One development is to find out that the strong 
tensor force is of intermediate range and we are able 
to express the tensor correlation in a reasonable shell 
model space pd| - [T3j . We name this method as Ten- 
sor Optimized Shell Model (TOSM). The other is the 
Unitary Correlation Operator Method (UCOM) to treat 
the short-range correlation caused by the short-range re- 
pulsion (l4Hla|. We shall combine these two methods, 
TOSM and UCOM, to describe nuclei using a bare NN 
interaction. In the TOSM part, the wave function is 
constructed in terms of the shell model basis states with 
full optimization of two particle-two hole (2p2h) states. 
There is no truncation of the particle states within the 
2p2h space in TOSM, where the spatial shrinkage of the 
particle states is essential to obtain convergence of the 
tensor contribution involving high momentum compo- 
nents 

EMU- 

This treatment of the bare tensor force in 
TOSM corresponds to the one-pair approximation corre- 
lated by the tensor force |19l - [2l| . In a few-body frame- 
work, the validity of TOSM was confirmed by taking only 



2 



the -D-wave component connected with the S'-wave state 
directly by the tensor force for the s-shell nuclei plj . 
Their calculation reproduces more than 90% of the ten- 
sor contribution and this result justifies the description 
of the tensor correlation in TOSM. The explicit inclu- 
sion of 2p2h states in the extended mean field model for 
heavy nuclei has been formulated by Ogawa et al. as an 
extended Brueckner Hartree-Fock theory [22j • 

So far, we obtained various successful results of 
TOSM for the investigation of the tensor correlations 
in He isotopes. In 4 He, we showed the selectivity of 
(Pi/2) 2 ( s i/2) -2 configuration in the 2p2h space with the 
pn pair induced by the tensor force and this correlation 
was recognized as the deuteron-like state [l2|. This spe- 
cific 2p2h excitation plays a decisive role to reproduce 
the spectra of neutron-rich He isotopes [l3|, [24| ■ It 
is found that the p 3 / 2 occupation of extra neutrons in- 
creases the tensor correlation of nuclei from 5 Hc to 8 He, 
while the P1/2 occupation of extra neutrons decrease the 
tensor correlation of those nuclei due to the Pauli block- 
ing between the specific 2p2h excitations by the tensor 
force in 4 He and the motions of the extra px/2-neutrons. 
This configuration dependence of the tensor correlation 
produces the right amount of splitting energy between 
the P3/2- and p^-dominant states in the He isotopes. 

In this paper, we use the TOSM+UCOM for the study 
of the Li isotopes from 5 Li to 9 Li, and discuss their struc- 
tures focusing on the role of the tensor force on the exci- 
tation energies and configurations of those nuclei. Before 
the present analysis, for neutron-rich 1 Li and 11 Li, we 
have performed the coupled two- and three-body cluster 
model analyses assuming the 9 Li core described usin g th e 
simplified TOSM wave function, respectively @, [25l. |26| . 
The cluster model for two nuclei precisely describes the 
neutron wave functions in a loosely bound state around 
9 Li, such as a neutron halo structure in n Li. Similar to 
the case of the He isotopes, we have shown that 9 Li has 
the amount of the selected excitation from the 0s orbit 
to the Opi/2 orbit by the tensor force, and the P1/2 occu- 
pation of last one and two neutrons in 10 Li and n Li are 
blocked by this tensor correlation in 9 Li. As a result, the 
Is occupation of last neutrons is relatively enhanced in 
two nuclei, and the virtual s-state in 10 Li, and the large 
s-wave mixing and a neutron halo formation in 11 Li are 
naturally explained. On the basis of these successful re- 
sults for two nuclei, in the present study, we perform 
a full microscopic analysis by using the TOSM+UCOM 
approach using the NN interaction for the Li isotopes. 
We investigate the systematic role of the tensor force in 
5_9 Li. It is interesting to see the effect of tensor force on 
the structures of 5 ' 6 Li [III, El]. In 6 Li, the isospin T = 
and T = 1 states coexist in the low excitation energy 
region. The influence of the tensor force on two kinds 
of states in 6 Li are important to understand the isospin 
dependence of 6 Li in relation with the tensor correlation. 
In heavier Li isotopes, the contribution of tensor force is 
investigated in the ground and excited states. We also 
see the high momentum behavior of the kinetic energy 



in each state in relation with the tensor correlation. The 
present subject of Li isotopes becomes a foundation of the 
nuclear structure study using TOSM+UCOM and also of 
the previous study of n Li and 10 Li structures [26j. 

In Sec. HU we explain the method of the 
TOSM+UCOM approach. In Sec. [Hi we show the 
level structures of Li isotopes and discuss their charac- 
teristics in relation with the tensor force. A summary is 
given in Sec. [TV] 



II. THEORETICAL FRAMEWORK 

A. Tensor-optimized shell model (TOSM) 

We explain the tensor optimized shell model (TOSM) 
for open shell nuclei. We begin with writing a many-body 
Hamiltonian for a A body system as, 

A A 

H = Y^T i -T c . ia .+J2v ij , (1) 

V V = vC+vT+vtf + v^. (2) 

Here, T, is the kinetic energy of each nucleon with T c m . 
being the center of mass kinetic energy. We take a 
bare interaction Vij such as AV8' Q consisting of cen- 
tral tensor vTa and spin-orbit vfj terms and vf/ mb 

t J L J * t J <>3 

the Coulomb term. Wc obtain the many-body wave func- 
tion '5 by solving the Schrodingcr equation H'fy = Ety . 
In our previous works for 4 He and a few body systems, 
we found that the tensor force can be described by tak- 
ing 2p2h excitations with lar ge m omentum components 
in the shell-model framework [12l. [2lj . 

We prepare first a standard shell-model state for an 
open shell Li nucleus with A nucleons in order to in- 
troduce the TOSM for open shell nuclei. The standard 
shell-model state is written as 

*s = J2 A ^\M 4 (0p) A -^ks) ■ (3) 

ks 

Here, the p shell is the valence shell and the index k$ is 
used to distinguish various shell-model components. The 
spirit of the TOSM is that the tensor force works strongly 
for two nucleons in the standard shell-model state and ex- 
cite two nucleons to various two particle states with high 
momentum components. Hence, we limit configurations 
up to two particle-two hole excitations from the standard 
shell-model state. 

We extend therefore first the standard shell-model 
state by allowing two particles in the s shell to excite 
into the p shell. We name these extended shell-model 
states as |0p0/i; feo) states, since no particle is excited into 
shell-model orbits outside of the sp shells. We write the 
extended shell-model states as 

\0p0h;k ) = \(0 S ) n °(0p) n *>;k ) (4) 
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with the constraints n s + n p = A and 2 < n s < 4. fco dis- 
tinguishes various sp shell configurations. We take now 
all the necessary high momentum components brought by 
the strong tensor force. These components are included 
by exciting two nuclcons from the sp shells to higher 
shells and therefore limit these configurations to 2p2h 
states from the standard shell-model state ^s- Hence, 
we have 



\2p2h;k 2 ) = |(0s) n "(0p) n » (higher) 2 ; k 2 ) 



(5) 



with the constraint n s + n p = A — 2 and 2 < n s < 4. 
We introduce quantum numbers k 2 to distinguish various 
2p2h states, which amount to a large number of configu- 
rations. Here, "highe.^ indicates higher shells outside of 
the sp shells, which arc treated as particle states. Hence, 
we use the notation of 2p2h states to specify that two 
particles are in higher shells outside of the sp shells. We 
allow also \p\h excitations for shell model consistency. 



|lpl/i;fci> = \(0s) n °(0p) n r (higher) 1 ;h) 



(6) 



with the constraint n s + n p = A — 1 and 2 < n s < 4. 
These lplh states are able to include high momentum 
components and improve the OpOh wave functions in the 
presence of the strong tensor correlation. Altogether, we 
write the TOSM wave function ^ for open shell nuclei as 



* = ^A^lOpQhiktti + ^AfiAlplhiki) 

ko ki 

+ J2^ 2 \2 P 2h;k 2 ). (7) 



Here, all the amplitudes {Ak Q , A^ , Ak 2 } are variational 
coefficients to be fixed by the energy minimization. 

We write now the details of the radial wave functions. 
The OpOh states are shell model states and expressed 
in terms of harmonic oscillator wave functions. Hence, 
the Os and Op shell-model wave functions are used to ex- 
press the radial wave functions, whose length parameters 
are taken independently as variational parameters. The 
lplh and 2p2h states involve particle states with high 
momentum components. The hole states correspond to 
the shell-model states in the sp shells. The particle wave 
functions have to contain high momentum components 
to express the specific characters of the tensor force with 
all the possible angular momenta until the total energy 
converges. 

For particle states, we employ the Gaussian wave func- 
tions to express single-particle states in higher orbits in 
order to describe high momentum properties due to the 
tensor force [29|, H3| • We prefer the Gaussian wave func- 
tions over the shell-model states in order to effecti vely 
include the necessary high momentum components [l2j. 
When we superpose a sufficient number of Gaussian wave 
functions, the radial components of the particle states 
can be fully expressed. Gaussian basis states should 
be orthogonalized to the hole states and among them- 
selves. This condition is imposed by using the Gram- 
Schmidt orthonormalization. In order to use the non- 
orthogonal Gaussian basis functions in the shell model 



framework, we construct the following orthonormalized 
single-particle basis function ij}™ using a linear combina- 
tion of Gaussian bases {4>a} with length parameter b aM . 



E 

jy=l 



CM = ><,u <t> a (r,b a , u ), 



<CIC<) 

for n 



1,. 



■,N n 



(8) 
(9) 



where N a is a number of basis functions for the orbit a, 
and v is an index to distinguish the bases with Gaussian 
length of b a>lJ . The explicit form of the Gaussian basis 
function is written as 



Ni{b a , v ) 



o /,-(2i+3) 

r(; + 3/2) 



(ii) 



where I and j are the orbital and total angular momenta 
of the basis states, respectively. The weight coefficients 
{d™ „} are determined to satisfy the overlap condition 
in Eq. ©. This is done by using the Gram-Schmidt 
orthonormalization. Following this procedure, we obtain 
the new single-particle basis states {f/'S} m Eq. © used 
in TOSM. The particle states in lplh and 2p2h states 
are prepared to specify the basis wave functions, whose 
amplitudes are determined by the variational principle. 

We note that we may use another method to obtain 
{ipa} by solving the eigenvalue problem of the norm ma- 
trix of the Gaussian basis set in Eq. (jTU)) with the di- 
mension N a . This method gives different coefficients 
Wa u) horn those of the Gram-Schmidt method. How- 
ever, these two methods of making the orthonormalized 
single-particle basis states provide equivalent variational 
solutions for the total TOSM wave function ^> in Eq. (|7|), 
because we start from the same Gaussian basis functions 
in Eq. (|10p with the same number N a of Gaussian basis 
states. Therefore, the TOSM solution does not depend 
on the construction of the orthonormal basis states when 
we start from the Gaussian basis functions. 

We construct Gaussian basis functions of particle 
states to be orthogonal to the occupied states. For the Is 
orbit in the particle states, we prepare an extended Is ba- 
sis function [TJ [l2j , which is orthogonal to the Os state 
and possesses a length parameter bi s ^ that can differ 
from bo s of the Os state. In the extended Is basis func- 
tions, the polynomial part is changed from the usual Is 
basis states to satisfy the conditions of the normalization 
and the orthogonality to the Os state [ll[ ■ For the lp or- 
bits, we take the same method as used for the Is case. In 
the numerical calculation, we prepare 10 Gaussian basis 
functions at most with various range parameters to get a 
convergence of the energy and Hamiltonian components. 

We note here that when we write probabilities and oc- 
cupation numbers of each orbit in various states in the 
numerical sections, those numbers are given by the sum- 
mation of all the orthogonal orbits with the same spin 
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having different radial behaviors due to the fact these 
wave functions are constructed by the orthonormaliza- 
tion. For hole states, we do not have to sum up their 
numbers because the hole states are described by the 
single harmonic oscillator basis states. 

We take care of the center-of-mass excitations by us- 
ing the well-tested method of introducing a Hamiltonian 
of center-of-mass motion in the many-body Hamiltonian 
known as the Lawson method [3l| . In the present study, 
we take the value of fouj for the center of mass motion 
as the averaged one used in the Os and Op orbits in the 
OpOh states with the weight of the occupation numbers in 
each orbit [l~3|. Adding this center of mass Hamiltonian 
as the Lagrange multiplier to the original Hamiltonian 
in Eq. (flj, we can effectively project out only the lowest 
harmonic oscilattor state for the center-of-mass motion. 

The variation of the energy expectation value with re- 
spect to the total wave function ^ in Eq. ([7]) is given 
by 



. 



which leads to the following equations: 
d{<S/\H - E\V) 



Sb a ,v 
d(^\H-E\^) 



OA 



= 
= 



ki 



for i = 0,1,2. 



(12) 

(13) 
(14) 



The total energy is represented by E. We solve two vari- 
ational equations in Eqs. (fT3j) and (fLlj) in the following 
steps. First, fixing the length parameters b a l , and the 
partial waves of the basis states up to L m ax, we solve the 
linear equation for {Ai ii } as an eigenvalue problem for H . 
We thereby obtain the eigenvalue E, which is functions 
of {b a ,v} and L max - Next, we try to adopt various sets of 
the length parameters {b a , v } and increase L max in order 
to find a better solution which minimizes the total en- 
ergy E. In TOSM, we can describe the spatial shrinkage 
of particle states with an appropriate radial form in the 
individual configuration, which is important to describe 
the tensor correlation [ll[ , as seen in the deuteron. 



B. Unitary Correlation Operator Method (UCOM) 

We briefl y expla in UCOM for the short-range central 
correlation p^ - tlq . in which the following unitary oper- 
ator C is introduced 



C = exp(-i^2gij) . 



(15) 



'<j 



We express the correlated wave function 'J in terms of a 
simple wave function $ as ^ = CQ. The transformed 
Schrodinger equation becomes HQ = EQ where the 
transformed Hamiltonian is given as H = C^HC. The 
operator C is in principle a many-body operator. In case 



TABLE I: Probabilities of various configurations in 4 He with 
TOSM+UCOM using the AV8' interaction, where the two 
subscripts 00 and 10 are the spin-isospin quantum numbers. 



(0s) 



(Os)oo 

10 2 (OPl/2)?0 



(o S )r 2 [(is 1/2 )(od 3 /2)]io 

(0 S )r 2 [(0p3/2)(0/ 5 /2)]lO 

(o s )r 2 [(op 1/2 )(op3/ 2 )]io 



remaining part 



84.14 
2.32 
2.20 
1.82 
1.21 



8.31 



of the short-range correlation, we are able to truncate 
the modified operators at the level of two-body opera- 
tors Q. 

Two-body Hermite operator g in Eq. (|15|) is defined as 



1 
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{p r s{r) + s(r)p r } 



(16) 



where the operator p r is the radial component of the 
relative momentum and is conjugate to the relative co- 
ordinate r. The function s(r) expresses the amount of 
the shift of the relative wave function at the relative co- 
ordinate r for every nucleon pair in nuclei. We use the 
TOSM basis states to describe <f>. 

In UCOM, the function s(r) is determined variation- 
ally to minimize the total energy of the system. We 
parametrize s(r) in the same manner as proposed by 
Feldmeier and Neff [3, EH for four channels of spin- 
isospin pair. The detailed forms of s(r) and their 
parametrization are explained in the previous paper fl3T ] . 
In the present analysis, we use these s(r) functions com- 
monly for every states of the Li isotopes. To simplify the 
numerical calculation, we adopt the ordinary UCOM for 
the central correlation part instead of the S'-UCOM in 
this analysis. 



III. NUMERICAL RESULTS 
A. 4 He 

We explain first the results of 4 He calculated by using 
TOSM+UCOM. We then discuss how the structures of 
5 ~ 9 Li behave by comparing with the properties of 4 He. 
We show results of 4 He using the AV8' interaction, which 
consists of central, LS and tensor terms. The AV8' inter- 
action is used in the calculation given by Kamada et al, 
where the Coulomb term is ignored ■ The total energy 
of 4 Hc is obtained as -22.30 MeV in TOSM+UCOM. 
The convergence of the solutions of 4 He has already been 
confirmed in the previous works [l2l . The energy 

components of He are discussed in the previous pa- 
per [l| 

The dominant configurations of 4 Hc in TOSM are 
listed in Table U with their probabilities. The specific 
2p2h states such as (0s) 1 " 2 (0p 1 / 2 )io show large probabil- 
ities and are essential to produce the tensor correlation 
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TABLE II: Occupation numbers in each orbit of 4 He using 
AV8' and Minnesota (MN) interactions. 



"He^) 


OSl/2 


Pl/2 


P3/2 


lSl/2 


0*3/2 


0*5/2 


AV8' 


3.77 


0.06 


0.04 


0.03 


0.04 


0.01 


MN 


3.94 


0.01 


0.03 


0.01 


0.004 


0.006 



TABLE III: Matter radii of Li isotopes in comparison with 
the experiments; aj3(§], b[37]]. Units are in fm. 





TOSM 


Experiment 


6 Li 


2.22 


2.35(3) a 2.44(7) b 


7 Li 


2.32 


2.35(3) a 


8 Li 


2.42 


2.38(2) a 2.50(6) b 


9 Li 


2.46 


2.32(2) a 2.44(6) b 



in 4 He. The selectivity of those 2p2h configurations can 
be understood from the coupling by the tensor opera- 
tor which changes both the orbital angular momentum 
and spin by two with the opposite direction of their z- 
components [U, 0). In addition, the 2p2h state with 
(0s) u0 2 [(ls 1 / 2 )(0c?3/2)]io is also favored by the tensor in- 
teraction. It is found that the dominant 2p2h states com- 
monly correspond to the excitations of a pn pair. Those 
features of the 2p2h excitations play an important role 
to determine the structures of 5 ~ 9 Li as will be discussed 
later. 

For comparison, we show the results with the Min- 
nesota (MN) NN interaction, which does not have the 
tensor force. The MN interaction consists of the central 
and LS parts with u parameter being 0.95 [HI, [33|. We 
use set III of the MN interaction for the LS force. TOSM 
has been confirmed to reproduce the binding energy and 
the radius of 4 He obtained by few-body stocastic varia- 
tional method [l3l. l3~4l |35| . It is obtained that the (0s) 4 
configuration is by far dominant in the 4 He wave function 
by 96.6%. Other 2p2h configurations have the probabil- 
ity less than a percent. This comparison makes the role of 
the tensor force clear on the energy spectra of the p-shell 
nuclei for the case of the AV8' interaction. 

We also list the occupation numbers of nucleon in 4 He 
using AV8' in Table ITT1 It is shown that the pi/2 orbit 
has the largest contribution among the particle states 
according to the large 2p2h mixing including the P1/2 
component. In the MN case, it is found that the com- 
ponent of the 0s orbit is larger than the AV8' case and 
the enhancement of the p 1 / 2 orbit is not obtained. These 
results mean that the tensor force brings the specific ex- 
citations from the s-shell to the p- and sd-shells in 4 He. 



B. Energy spectra of the Li isotopes 

We show the calculated results of the Li isotopes using 
TOSM+UCOM with the AV8' interaction, where i raax 
is taken as 10 to get sufficient convergence. We show the 
excitation spectra of 5 ~ 9 Li in Fig. [TJ We sec quite a good 
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FIG. 1: Excitation energies of Li isotopes with 
TOSM+UCOM using AV8'. 
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FIG. 2: Excitation energies of He isotopes using AV8' with 
TOSM+UCOM. 



correspondence to the experimental spectra. Similar to 
this result, we obtain good results for the excitation en- 
ergies in the He isotopes as shown in Fig. [2] the details of 
which were discussed in Ref. [l3| . In Fig. [TJ the resulting 
level spacing of the Li isotopes in TOSM+UCOM is good, 
but slightly more compact than the experimental spectra. 
For example, in 9 Li, the small energy difference between 
the lowest 3/2" and 1/2" states in TOSM+UCOM in 
comparison with the experiment. These characteristics 
are commonly obtained in the GFMC calculation (38| . 
We will discuss in detail all the level structures of the 
Li isotopes in the next sub-section. The additional gen- 
uine three-body interaction can be one of the components 
to reproduce the experimental situation jij. As for to- 
tal binding energies, our results underestimate them in 
He and Li isotopes and its amount becomes larger for 
neutron-rich side. 

To see the interaction dependence, the energy spectra 
using the MN interaction are shown in Fig. [3] In the 
MN case, the excitation energies arc split in two regions 
in comparison with the experiments. One of the reasons 
is the too large splitting energy between p 3 / 2 and P1/2 
components mainly given by the LS interaction. When 
we adjust the MN interaction such as changing the u pa- 
rameter and the strength of the LS interaction, to fit the 
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FIG. 3: Excitation energies of the Li isotopes using the Min- 
nesota interaction (MN) with TOSM. 




4 He 6 He "He 



FIG. 4: Matter radius (R m ) of the He and Li isotopes in 
TOSM+UCOM with solid circles. Open circles are the ex- 
periments; [39I-I4H for He isotopes and [36l.l37| for Li isotopes. 



LS splitting energy of 5 Li, it is still difficult to reproduce 
the whole trend of the energies of the ground and excited 
states of Li isotopes consistently. A large difference be- 
tween the AV8' and MN interactions is the tensor force. 
The comparison of the energy spectra of the two interac- 
tions indicates the effect of the tensor force, which has a 
large impact on the excitation energy spectra of the He 
and Li isotopes. 

We show total energies of the Li isotopes as shown in 
Fig. [5] The experimental total energies become lower 
than that of 4 He as the neutron number increases. On 
the other hand, the calculated results go up in energy 
with the neutron number. Hence, we lack some attrac- 
tions, which increases with the mass number. This luck of 
mass dependent attraction was already seen in the case 
of the He isotopes. In the previous analysis of the He 
isotopes [13j . we have already discussed the reasons of 
undcrbinding, which can contribute to the bulk part of 
the binding energies; 

1) higher configurations beyond the 2p2h states in 
TOSM, such as the ApAh states to include the two-kinds 
of the 2p2h excitations with isospin T = and T = 1 
pairs simultaneously. 

2) the genuine three-body interaction, and 
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FIG. 5: Total energies of Li isotopes using AV8' measured 
from that of 4 He. The experimental energy spectra are shown 
in the left hand side of the corresponding theoretical spectra. 



3) the improvement of the correlation function form of 
R+(r) in UCOM. 

We have estimated these three effects on the binding 
energies of He isotopes such as 6 He and 8 He and found 
that we can obtain sufficient binding energies for He iso- 
topes, which are comparable with the experimental val- 
ues [13| . In the present study, from the results of energy 
spectra in Figs. [T]and[2l we focus on the discussion of the 
structure differences between individual energy levels of 
the Li isotopes. We show also the matter radii of 6 > 7 ' 8 > 9 Li 
in TOSM+UCOM in Table HIl which agree with experi- 
ments. We include the results of He isotopes in addition 
to the Li case in Fig. |4]and find the whole trend on the 
matter radii observed in the He and Li isotopes is very 
good. 



C. Energy components and configurations 



'Li 



We discuss in detail the structures of each energy level 
of the Li isotopes obtained by using the TOSM+UCOM 
framework with the AV8' NN interaction. In Fig. [TJ 
we show the energy spectrum of 5 Li as compared with 
experiment, where the energy difference of the spin-orbit 
partner states is obtained as AE = 1.7 MeV. It is noted 
that in the present calculation, the continuum effect of 
the last unbound proton is not included because of the 
bound state approximation. It is naively expected that 
the inclusion of the continuum state of a last proton can 
reduce the splitting energy of 5 Li because of the spatial 
extension of the proton wave function in the asymptotic 
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TABLE IV: Various energy components in 5 Li measured from 
those of the 4 He ground state. Units are in MeV. 

J 71 Kinetic Central Tensor LS 
3/2- 10.96 -2.00 -3.07 -1.74 
1/2- 5.90 0.30 -0.41 0.03 

TABLE V: Various energy components in 6 Li measured from 
those of the 4 He ground state. Units are in MeV. 

,F(T) Kinetic Central Tensor LS 



It 


(0) 


28.29 


-9.76 


-11.81 


-2.28 


3+ 


(0) 


27.57 


-8.75 


-8.75 


-3.35 


2+ 


(0) 


25.90 


-8.08 


-8.69 


-1.29 


1 2 + 


(0) 


23.96 


-6.63 


-6.44 


-1.32 


0+ 


(1) 


25.72 


-13.35 


-4.12 


-2.29 


2+ 


(1) 


30.50 


-12.00 


-6.64 


-3.07 



region. There is also a discussion of the effect of the 
genuine three-body interaction on the splitting energy in 
the 4 He+n scattering [Hj]. 

In Table IIV1 we compare various energy components 
in the 3/2" and 1/2 - states of 5 Li measured from those 
of 4 He. We discuss the double roles of the tensor force on 
these spin-orbit partner states. From Table llVl it is found 
that the tensor energy of the 3/2 _ state is attractive, 
while it is very weak for l/2 — . The larger contribution 
of the tensor force in 3/2 _ results in the enhancement 
of the kinetic energy because of the involvement of high 
momentum components brought by the tensor force. For 
1/2 - , on the other hand, the energy gain from the tensor 
force is small and the enhancement of the kinetic energy 
is small. These differences are related to the large mixing 
of the pi/2 component in 4 He as shown in Tablcfl] In 5 Li, 
when the last proton occupies the p^/2 orbit, the ten- 
sor correlation in 4 Hc is not much reduced. On top the 
P3/2 state has an additional tensor contribution because 
of the tensor correlation between the p 3 / 2 proton and the 
neutrons in the s-shell. In case of the P1/2 occupation, 
this proton blocks the spatially compact component of 
the p\/2 proton in the 4 He part because of the small de- 
generacy of the pi/2 orbit. This dynamics produces the 
Pauli-blocking and reduces the total binding energy of 
5 Li. As a result, the last proton occupied in the P1/2 
orbit should be orthogonal to the excited P1/2 orbit in 
4 He and the tensor force does not gain the energy in 
5 Li(l/2 _ ). Those different couplings of the tensor force 
are essential to explain the energy components in two 
states of 5 Li and results in the LS splitting energy as a 
net value dHHSl 

2. 6 Li 

In Fig. [TJ we show the energy spectrum of 6 Li as com- 
pared with experiment. We see the ground state has the 



spin-parity of J x = 1 + , while states with T = 1 are closer 
to the 1 + state as compared to experiment. In order to 
understand the reason of the calculated results, we show 
energy components in Table [V] measured from the energy 
components in 4 He. It is found that the 1 + ground state 
mostly exhausts the tensor energy. This is because of 
the deuteron-like correlations of the valence pn-pair in 
addition to that in the 4 He part. Corresponding to this 
fact, the kinetic energy of the if state shows the largest 
value among the isospin T = state of 6 Li. The 2 + and 
3 + states with T = almost show the similar energies 
for the central and tensor energy components. For the 
two isospin T = 1 states, corresponding to the isobaric 
analog states of 6 He, the tensor energies are mostly much 
smaller than the T = states of 6 Li. This is naturally 
understood from the viewpoint of the isospin dependence 
of the tensor force. In the T = 1 states of 6 Li, the 4 He 
part contains the pn pair excitation with T = from 
the Os-shell by the tensor force, but the outer pn pair in 
the Op-shell with a T = 1 state shows a small amount of 
the excitation due to the weak tensor force in the T = 1 
channel. The 2 + state with T = 1 gains the larger tensor 
energy than those of + state and show the larger kinetic 
energy. 

For energy levels of 6 Li in comparison with experiment, 
the TOSM results give a small energy difference between 
T = and T = 1 states as shown in Fig. Q] This can be 
related with the clustering effect of a+d structure of 6 Li 
with T = 0. The TOSM wave function is expanded in 
terms of the shell model basis states and hence, is naively 
difficult to describe the tail component of the spatially 
extended clustering states, such as the 3a state in the 
excited states of 12 C [H. Similarly, 6 Li with the T = 
states is considered to need the clustering component of 
a+d in addition to the shell model like one [13, HH. The 
lacking of the cluster components loses some amount of 
the binding energy of 6 Li in TOSM. On the other hand 
the T = 1 states do not corresponds to the a+d clustering 
state, so that TOSM can describe those T = 1 states, 
similar to the 6 He case. This structure difference between 
T = and T = 1 states can explain the small energy 
difference between those states of 6 Li in TOSM. When 
we increase the strength of tensor force artificially by 
about 5% in order to see the effect of tensor force on 
the T = and T = 1 states in 6 Li, the energy gains of 
the 1 + (T = 0) state is 3.1 MeV, which is larger than 2.7 
MeV of the + (T = 1) state. This fact implies that the 
T = states in 6 Li contains the tensor correlation much 
more due to the presence of the outer pn pair with the 
T = channel. The full inclusion of the a+d component 
in addition to the TOSM basis states would be one of 
the ways to explain the experimental energy difference 
between the T = and T = 1 states of 6 Li. 

We list the dominant configurations of 6 Li in Table IVT1 
For comparison of the if and \\ states, it is found that 
the if state does not have the large component of the 
(Opi/2) 2 configuration of a pn pair, while the \\ state 
has this configuration by 23%. This difference in the 
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TABLE VI: Dominant configurations of 6 Li( J*"(T)) with their 
squared amplitudes (A'/.) 2 using AV8' interaction. 



tf(0) 


l + (0) 


(0s) 4 (0p 1/2 )(0p 3/2 ) 


0.43 


(0s 


) 4 (0p 3 /2) 2 


0.30 


(0s) 4 (0p 3 / 2 ) 2 


0.38 


(0s 


) 4 (0p 1/2 )(0p 3/2 ) 


0.29 






(0s 


) 4 (0pi/ 2 ) 2 


0.23 



2+(0) 


3+(0) 


(0s) 4 (0p 1/2 )(0p 3/2 ) 


0.82 


(0s) 4 (0p 3/2 ) 2 0.82 



TABLE VIII: Occupation numbers in each orbit of Li using 
AV8' interaction. 





0sl/2 


0pi/ 2 


0P3/ 2 


lSl/ 2 


^3/ 2 


d 5 /2 


1+0 


3.74 


0.53 


1.42 


0.04 


0.05 


0.04 


i + ,o 


3.74 


0.89 


1.08 


0.05 


0.05 


0.04 


2+,0 


3.75 


0.97 


0.98 


0.04 


0.05 


0.04 


3+,0 


3.74 


0.03 


1.94 


0.05 


0.05 


0.03 


o+,i 

2+,l 


3.73 
3.74 


0.28 
0.15 


1.70 
1.83 


0.04 
0.04 


0.05 
0.05 


0.03 
0.03 



0+(l) 


2+(l) 


(0s) 4 (0p 3/2 ) 2 0.72 


(0s 


) 4 (0p 3/2 ) 2 


0.74 


(0s) 4 (0p 1/2 ) 2 0.11 


(0s 


) 4 (0p 1/2 )(0p 3/2 ) 


0.11 



TABLE VII: Dominant configurations of 6 Li(J* (T)) using 
Minnesota interaction. 



l+(0) 


i 2 + (o) 


(0s) 4 (0p 3/2 ) 2 


0.56 


(0s 


) 4 (0pi/ 2 )(0p 3 /2) 


0.43 


(0s) 4 (0p 1/2 )(0p 3/2 ) 


0.29 


(0s 


) 4 (0p 3/2 ) 2 


0.19 






(0s 


) 4 (0p 1/2 ) 2 


0.11 



2+(0) 


3+(0) 


(0s) 4 (0p 1/2 )(0p 3 / 2 ) 


0.83 


(0s) 4 (0p 3/2 ) 2 0.88 



0+(l) 


2+(l) 


(0s) 4 (0p 3/2 ) 2 0.82 


(0s 


) 4 (0p 3/2 ) 2 


0.77 


(0s) 4 (0p 1/2 ) 2 0.04 


(0s 


) 4 (0pi/ 2 )(0p 3/2 ) 


0.09 



configurations explains the tensor energies in two states 
shown in Table [V] Owing to the Pauli blocking between 
the pn-pair in the Op-shell and another pair excited from 
the Os-shell, the 1+ suffers the large blocking effect due to 
the presence of the (Opi/2) 2 configuration, which reduces 
the tensor contribution in this state. The same state- 
dependence of the tensor correlation is mentioned in 5 Li. 

We also compare the dominant configurations of 6 Li us- 
ing AV8' and MN interactions. The MN case is shown in 
Table \VU\ In the 1+ ground state, the difference in the re- 
sults of two interactions is clearly seen; the (0pi/ 2 )(0p 3 / 2 ) 
configuration in the Op-shell of 6 Li is enhanced in the 
AV8' case rather than the MN one. This configuration is 
a characteristic feature of the LS coupling scheme. We 
discuss the relation between this configuration and the 
tensor force. The (0pi^)(0p3/2) configuration of pn in 
6 Li involves the 89% of the spin triplet (S = 1) com- 
ponent, while the (OP3/2) 2 configuration involves 44%. 
The S = 1 component in the wave function is impor- 
tant to activate the tensor force and is related to the 
deuteron-like component in the pn pair. It is found that 
the (Opi/2) (OP3/2) configuration in 6 Li is favored to in- 



TABLE IX: Occupation numbers in each orbit of Li using 
MN interaction. 



J\T 


0si/ 2 


0pi/ 2 


0p3/ 2 


lsi/ 2 


^3/ 2 


d 5 /2 


l+,0 


3.85 


0.37 


1.60 


0.03 


0.03 


0.04 


1+0 


3.83 


0.76 


1.08 


0.05 


0.04 


0.06 


2+,0 


3.85 


0.94 


0.99 


0.03 


0.04 


0.04 


3+,0 


3.89 


0.02 


1.92 


0.04 


0.01 


0.04 


0+1 

2+,l 


3.85 
3.85 


0.12 
0.12 


1.85 
1.84 


0.03 
0.04 


0.02 
0.02 


0.05 
0.04 



crease the tensor energy in the wave function. Contrast- 
ingly, in the MN case, there is no explicit tensor force 
and the (0p 3 / 2 ) 2 configuration becomes dominant, which 
can be understood naively from the jj-coupling scheme. 
It is experimentally interesting to observe the amount of 
the mixings of P1/2 and P3/2 orbits in 6 Li. 

We list the occupation numbers of 6 Li up to the sd- 
shell in Table ETTTl using AV8'. It is found that the oc- 
cupation numbers in the Os-shell and the sd-shell show 
almost common values among the 6 Li states, and the 
Op-shell contributions depend on the states. For compar- 
ison, we also list the MN case in Table IIXI It is found 
that the 0s components in each state of 6 Li using AV8' 
is larger than the MN case. This is similar to the 4 He 
results shown in Table [TTJ Those results mean that the 
tensor force brings the excitation of the nucleon from 
the Os-orbit to higher orbits such as pi/ 2 in 6 Li. This 
tendency can be commonly seen in every state of 6 Li. 
Among the sd-shell components of 6 Li, it is found that 
the (I3/2 component is the largest one in most of the 
states using AV8'. This aspect cannot be seen in the 
MN case shown in Table IIXI , which shows the normal 
order of the jj coupling scheme. This difference comes 
from the tensor force, which selectively excites the nu- 
cleon in the 0p 3 / 2 orbit in the hole state to the d 3 / 2 orbit 
in the particle state, both of which has the same j quan- 
tum number. This reason is schematically explained in 
Fig. |51 in which the direction of the coupled orbital an- 
gular momentum and that of the coupled spin should be 
opposite in the tensor operator. Considering the transi- 
tion from the (OP3/2) 2 configuration by the tensor force, 
one of the favored components is naively estimated to be 
the (0d 3 / 2 ) 2 one with the distribution of Zi=Z 2 =2, where 
there are some fragment of the orbital angular momen- 
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FIG. 6: Schematic figure to explain the selectivity of the ten- 
sor coupling of two nucleons from the (p3/ 2 ) 2 configuration to 
the (d 3 / 2 ) 2 one. 



TABLE X: Various energy components in 7 Li measured from 
those of the 4 He ground state. Units are in MeV. 





Kinetic Central 


Tensor 


LS 


1/2- 


53.59 


-28.75 


-16.45 


-2.56 


1/27/ 


46.85 


-18.76 


-13.48 


-2.44 


3/2- 


53.77 


-29.46 


-15.62 


-3.24 


3/27 


46.95 


-18.15 


-14.66 


-2.26 


5/2- 


50.63 


-22.95 


-14.74 


-2.40 


5/2- 


49.35 


-19.56 


-14.10 


-4.54 


5/23 


46.45 


-17.58 


-13.02 


-2.64 


7/2- 


52.74 


-24.96 


-13.99 


-4.11 


7/2 2 - 


47.54 


-18.02 


-12.82 


-3.76 



turn distribution of two nucleons in the particle states. 

The same selectivity of the tensor coupling can be seen 
in the case of 0si/ 2 -0pi/ 2 combination in 4 He. These spe- 
cific couplings can be seen from the view point of the pi- 
onic correlation in the particle-hole representation. The 
pion has the spin-parity of - . Hence, hole configura- 
tions with positive or negative parity are favorably ex- 
cited to particle configurations with opposite parity with 
the same total spin due to the pion exchange. The selec- 
tivity of the tensor force can be seen as the selectivity of 
the pion exchange [f| [Tt], HH . It is extremely interesting 
to observe experimentally the relatively large mixing of 
d 3 /2 orbit among the sd-shell in 6 Li. 



TABLE XI: Dominant configurations of 7 Li(J 7r (T)) with their 
squared amplitudes (A'j,) 2 ■ Two subscripts are the spin- 
isospin quantum numbers. 



1/27 


l/2 2 


(08 


) 4 (0p 3 /2)gi(0p 1/2 ) 


0.50 


(Os) 4 (Op3/ 2 )io(OPi/ 2 ) 


0.47 


(0s 


) 4 (0p 3 / 2 ) 3 


0.19 


(0 S ) 4 (0p 3/2 ) 3 


0.17 


(0s 


) 4 (OP3/2)io(OPi/2) 


0.07 


(0s) 4 (0p3/ 2 )oi(0pi/ 2 ) 


0.09 



3/2- 


3/27 


(0s 


) 4 (0p 3 / 2 ) 3 


0.48 


(0s) 4 (0p3/ 2 )!i(0p 1/2 ) 


0.55 


(0s 


) 4 (0p 3 /2)(0pi /2 )g 1 


0.15 


(Os) 4 (Op 3/2 )io(Opi/ 2 ) 


0.13 


(0s 


) 4 (0p3/ 2 )ii(0p 1/2 ) 


0.10 


(0s) 4 (0p3/ 2 )(0p 1/2 )g 1 


0.10 



3. 7 < 8 ' 9 £i 

For 7 Li, various energy components are shown in Ta- 
ble [X] measured from the 4 He components. From the 
comparison of the components in each state, it is found 
that the 1/2 7/ and 3/2 ~ states show the similar structure 
having the large tensor contributions and also the larger 
kinetic energy than those of other states. 

For comparison between the dominant configurations 
of 3/27/ and 3/2"/ states from Table IXll the 3/27/ state 
does not have the 0pi/ 2 orbit dominantly and this prop- 
erty increases the tensor contribution via the excitation of 
pn-pair from the Os-shell. On the other hand, 3/2,7 state 
have Opi/2 orbit in the most dominant configuration and 
this configuration blocks the excitation of pn-pair from 
the Os-shell and cannot increase the tensor contribution 
in comparison with the 3/27/ state. 

For comparison between the 1/27/ and 1/27/ states, the 
1/27/ state dominantly has the nn pair of the OP3/2 or- 
bit with the T — 1 component and the 1/277 state domi- 
nantly has the pn-pair of the 0p 3 / 2 orbit with T = 0. The 
latter configuration in the 1/277 state blocks the excita- 
tion of pn pair from the Os-shell to 0^3/2 orbit with some 
amount, which reduces the tensor contribution. These 
configurations in two l/2~ states determine the different 
tensor contributions in two states. From those analyses, 



5/27^ 5/2^ 

(0 S ) 4 (0p 3/2 )ii(0pi /2 ) 0.64 (0s) 4 (0p 3 / 2 ) 3 0.61 
(0s) 4 (0p 3 / 2 )(0p 1/2 ) 10 0.09 (0 S ) 4 (Op 3/2 )io(Opi/ 2 ) 0.12 
(Os) 4 (Op 3/2 )io(Opi/ 2 ) 0.06 (0 S ) 4 (0p 3 /2)(0pi/ 2 )?o 0.05 



5/23 

(Os) 4 (Op 3 / 2 )io(Opi/ 2 ) 0.60 
(0s) 4 (0p 3/2 ) 3 0.11 
(0s) 4 (0p 3/2 )ii(0p 1/2 ) 0.06 



7/2/ 


7/2 2 


(0s 


) 4 (0p 3/2 ) 3 


0.53 


(06 


) 4 (op 3 / 2 )!o(opi /2 ) 0.53 


(0s 


) 4 (Op 3 / 2 )io(Opi/ 2 ) 


0.29 


(06 


) 4 (0p 3/2 ) 3 0.29 



it is found that the nucleons occupied in Op-shell in the 
lowest configuration play important role to explain the 
tensor contribution in each state of 7 Li. 

For 8 Li, various energy components are shown in Table 
IXHI measured from the 4 He components. The dominant 
configuration are listed in Table IXITT1 In 8 Li, the 2 + state 
is the ground state and has the largest tensor contribu- 
tion in this nucleus. This property can be understood 
from the tensor correlation; The 2 + state dominantly has 
the (OP3/2) 4 configuration instead of using the Opi/2 or- 
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TABLE XII: Various energy components in Li measured TABLE XIV: Various energy components in Li measured 



from those of the He ground state. Units are in MeV. 



from those of the He ground state. Units are in MeV. 





Kinetic Central 


Tensor 


LS 


0+ 


62.47 


-34.15 


-14.30 


-1.07 


it 


64.04 


-34.21 


-14.71 


-3.51 


' 


62.35 


-33.64 


-13.11 


-2.18 




61.69 


-31.86 


-13.17 


-1.89 




59.58 


-29.30 


-11.59 


-2.88 


2r 


72.08 


-38.69 


-19.30 


-4.69 


2* 


69.67 


-36.53 


-16.41 


-3.10 


o+ 
->, 


69.23 


-37.16 


-14.08 


-3.71 


3+ 


70.02 


-36.65 


-15.68 


-4.86 


4+ 


62.91 


-31.11 


-12.43 


-3.32 



,r 


Kinetic Central 


Tensor 


LS 


1/2- 


90.78 


-52.00 


-20.08 


-4.52 


3/2- 


91.49 


-53.35 


-18.62 


-5.72 


3/2- 


89.76 


-50.40 


-19.24 


-3.15 


5/2- 


90.66 


-51.47 


-18.42 


-4.20 


7/2- 


90.05 


-49.00 


-17.68 


-5.19 



TABLE XIII: Dominant configurations of 8 Li(J 7r (T)) with 
their squared amplitudes (A'l) 2 . Two subscripts are the spin- 
isospin quantum numbers. 



0+ 


4+ 


(Os) 4 (0 P 3/2)oi(Opi/ 2 )oi 


0.61 


(0s) 4 (0p 3/2 ) 4 


0.81 


(0s) 4 (0p 3/2 ) 3 (0p 1/2 ) 


0.20 







TABLE XV: Dominant configurations of 9 Li(J' r (T)) with 
their squared amplitudes (A'l) 2 . Two subscripts are the spin- 
isospin quantum numbers. 

1/2: 

(0s) 4 (0p 3/2 )^(0p 1/2 ) 0.67 

(Os) 4 (Qp3/ 2 )oi(Opi/ 2 ) 3 an 



3/gr 

0.46 
0.19 
0.07 



(0 S ) 4 (0p 3/2 ) 5 

(0s) 4 (0p 3 / 2 )^ /2 , 1/2 (0p 1/2 )g 1 

(0s) 4 (0p3/2)^/ 2 ,3/ 2 (0pi/ 2 )gl 



(0 S ) 4 (0p 3/2 ) 3 /2! 3 /2 (0p 1/2 ) 0.32 
(0s) 4 (0p 3/2 ) 4 0.25 
(0s) 4 (0p 3/2 )^ /2 , 1/2 (0p 1/2 ) 0.14 



^2 

(0 S ) 4 (0p 3/2 )^ /2 , 1/2 (0p 1/2 ) 0.50 (0 S ) 4 (0p3 /2 )! /2 , 3/2 (0p 1/2 )g 1 



3/2" 



(0s) 4 (0p 3/2 ) 3 /2 , 3/2 (0p 1/2 
(0s) 4 (0p3/ 2 )io(0pi/ 2 )oi 



°- u (0 S ) 4 (0p3 /2 ) 4 !(0p 1/2 ) 

j^f (0^) 4 (QP3/ 2 )| /2 ,3/ 2 (OPl/ 2 )fo 



0.38 
0.27 
0.08 



(0s) 4 (0p 3/2 )g 1 (0p 1/2 ) 10 0.27 

(0s) 4 (0p 3/2 )^ /2 , 1/2 (0p 1/2 ) 0.23 

(0 S ) 4 (0p3 /2 )^ /2 , 3/2 (0p 1/2 ) 0.14 

(0s) 4 (0p 3/2 ) 10 (0p 1/2 )g 1 0.13 



1 4 + 



(0s) 4 (0p 3/2 ) 4 0.27 
(0s) 4 (0p 3/2 )i 1 (0p 1/2 ) 10 0.23 
(0s) 4 (0p3 /2 )^ /2 , 3 / 2 (0pi/ 2 ) 0.16 
(0s) 4 (0p 3/2 ) 10 (0p 1/2 )g 1 0.06 



5/2 _ 


7/2~ 


(0s) 4 (0p 3/2 )li(0p 1/2 ) 0.57 
(0s) 4 (0p 3/2 )i /2 , 1/2 (0p 1/2 )g 1 0.13 


(0s) 4 (0p 3/2 )!i(0p 1/2 ) 0.80 



2+ 



(0 S ) 4 (0p 3/2 ^ 4 



/2)- 0.41 
(0s) 4 (0p3 /2 )3> 2 ,3 /2 (0p 1/2 ) 0.15 



(0s) 4 (0p 3/2 )li(0p 1/2 )g 1 0.14 



(0s) 4 (0p 3 / 2 )i /2 ,3/ 2 (0pi/ 2 ) 
(0s) 4 (0p 3/2 )3 /2il/2 (0p 1/2 ) 
(0s) 4 (0p3/ 2 )li(0pi/2 " 



\ 2 
/2J10 



0.35 
0.28 
0.08 



(0 S ) 4 (0p3/2)^ /2 , 1/2 (0p 1/2 ) 0.39 
(0s) 4 (0p 3/2 ) 3 /2! 3 /2 (0p 1/2 ) 0.20 
(0s) 4 (0p 3/2 ) 4 0.16 



(0 S ) 4 (0p 3/2 ) 4 
(0s) 4 (0p 3/2 )| (0p 
(0 S ) 4 (0p 3/2 ) 3 /2 , 1/2 (0p 1/2 ) 0.12 



01/2)01 



0.53 
0.15 



bit. This allows the selected excitation from the Os-orbit 
to the Opi/2 one in the 4 He part and increase the tensor 
correlation in this state relatively in comparison with the 
other states of 8 Li. The 3 + state also show the similar 
result of the larger tensor contribution. 

For 9 Li, various energy components are shown in Ta- 
ble [XTV] and the dominant configuration are listed in Ta- 



ble IXV1 This nucleus is important in relation with the 
structures of more neutron-rich systems of 10 Li and 11 Li 
and the breaking of the Op-shell magicity in that region. 
From Table IXVl in the ground 3/27" state, the nucleons 
in the Op-shell are dominated by the (0p 3 / 2 ) 5 configura- 
tion which is a sub-closed configuration for neutron part. 
Two neutrons among them can be excited into Opi/2 orbit 
with about 26%, which is the neutron pairing correlation 
in the Op-shell. This situation of the 3/27/ state can al- 
low the excitation of pn-pair from Os-shell, which emerges 
the tensor correlation in 9 Li. In our previous work, we 
have used this idea to explain the breaking of the neu- 
tron magicity in N = 8 in n Li and also the inversion 
phenomena of p-sd shells in 10 Li. In that work, the 2p2h 
excitations from the neutron sub-closed configuration of 
9 Li is taken into account as the pn and nn pair excita- 
tions. It is also found that the main component of the 
excited 3/27 state can correspond to the pairing excited 
one with respect to the ground 3/27" state. 

Seeing the whole structures from 6 Li to 9 Li, it is found 
that their ground states possesses the largest tensor en- 
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ergies, in which the configurations of Op-shell are con- 
structed to activate the two-kinds of excitations by the 
tensor force; one is from Os-orbit to the P\/2 orbit and 
the other is from Op-shell to higher shell such as sd- 
shell. In their ground states, the 0p\/2 component in 
their most dominant configurations is included only in 
6 Li. Other heavier systems contain the {0p3/2) A ~ 4 con- 
figuration largely, which is the jj-like one, to increase 
the tensor correlation from the Os orbit to the Opi/2 or- 
bit. From this result, the 6 Li nucleus is a specific one and 
this is considered to be related to the a+d component in 
this nucleus, as was mentioned. 

IV. SUMMARY 

We have developed a method to describe nuclei with 
bare NN interaction on the basis of the tensor opti- 
mized shell model with the unitary correlation opera- 
tor method, TOSM+UCOM. We have treated the tensor 
force in terms of TOSM, in which 2p2h states are fully 
optimized to describe the deuteron-like tensor correla- 
tion. The short-range repulsion in the NN interaction is 
treated by using the central correlation part of UCOM. 
We have shown the reliability of TOSM+UCOM using 
the AV8' interaction to investigate the structures of the 
Li isotopes. It is found that the excitation energy spectra 
are found quite consistent with the experimental spectra. 
When we employ the effective Minnesota interaction con- 
sisting of only the central and LS parts, the excitation 
energy spectra show quite a large deviation from experi- 
ment. 

It has been found that 4 He contains a relatively large 
amount of the pn pair in the P1/2 orbit due to the tensor 
force. This characteristics of the tensor force produces 
the state-dependence in the Li isotopes. In 5 Li, the 3/2~ 
state gains more tensor energy than the l/2~ case. The 
enhancement of the kinetic energy is also observed in the 
3/2~ state of 5 Li, which is brought by the tensor force. 
As a result, tensor force dynamically produces the state 
dependence in 5 Li and contributes to the LS splitting en- 



ergy in 5 Li. In 6 Li, the tensor force also makes the large 
mixing of the spin-triplet component in the configura- 
tion of the last pn pair in the ground state, which shows 
the LS coupling scheme. This component increases the 
tensor energy and is related to the a+d clustering con- 
figuration in 6 Li. This characteristics cannot be seen 
in the Minnesota interaction without the tensor force, 
which shows rather the conventional jj-coupling scheme. 
In 7 ' 8 ' 9 Li, the jj-like configurations, which include more 
0p3/2 states than the 0p!/ 2 states, can gain the tensor en- 
ergy, because of the allowance of the pn-pair excitation 
from the Os-orbit to Opi/2 orbit in the 4 He core part. As 
a conclusion of the roles of the tensor force in the Li iso- 
topes, the tensor energy depends on the configurations of 
nucleons occupied mainly in Op-shell in each state. This 
property is related with the amount of the excitation of 
pn pair from Os-shell by the tensor force. 

The amount of the high momentum component also 
depends on the tensor energy in each state. Observation 
of the high momentum component experimentally in fi- 
nite nuclei is desired in order to confirm the existence of 
the strong tensor correlation 0, . 

The tensor force would play important roles on various 
physical quantities such as charge radii, quadrupole and 
magnetic moments (47j in addition to the energy spectra. 
In our previous study of 11 Li, we discuss those quantities 
of 11 Li in detail [25|, [2(|. It is interesting to discuss the 
effect of the tensor force on these quantities in other Li 
isotopes using the TOSM+UCOM method. 
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